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ABSTRACT: One of the devastating effects of the earthquake is on liquefaction phenomenon site, which is one of the most important and 
most complex topics in seismic geotechnical engineering. Liquefaction is a phenomenon that occurs in loose, saturated sediments without 
stickiness in undrained conditions under the influence of waves caused by earthquakes or heavy static load. In order to cope with devastating 
effects of this phenomenon, it is essential to identify areas prone to liquefaction. This identification can be conducted as microzoning, in 
which risk capacity is determined in different areas. Thus, in this paper, the zoning map of Babol liquefaction risk will be provided. In this 
regard, a study was conducted on the soils in Babol and after examining different areas of the city, laboratory results and field studies of 
more than 50 boreholes in different areas with a depth of 20 m were analyzed for finding liquefaction and non-liquefaction segments. In this 
study, different approaches were used including Seed, Iwasaki, Haeri and Yasrebi, Chin & Zhang and Sewmez & Gocojlou procedures and 
finally, a computer program was written for examining and providing microzoning map of Babol liquefaction risk. 
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1. INTRODUCTION 

Since the dawn of mankind and human settlement under the roof, 
earthquake was always one of the risks that have killed the human 
from time to time. In recent years, discussion about the effects of 
earthquakes on soil has attracted many soil engineers and 
researchers. Among the phenomena related to seismic geotechnical 
risks, soil liquefaction is of great importance(Gill and Berg 1967). 

In this phenomenon, the saturated soil will be influenced by light 
stresses caused by earthquake and pore water pressure increases. 
Now if the speed of applying these loads is low so that the soil has 
an opportunity for drainage, pore water pressure will be amortized 
and liquefaction will occur. But if the speed of applying these loads 
is high, there is no opportunity for drainage, and the soil mass will 
tend to change volume and pore water pressure will increase so that 
it will be equal to the total stress applied to the soil mass(Maurer, 
Green et al. 2015). Therefore, the sand particles will be floating in 
the water and effective stress will be equal to zero. In this case, soil 
bulk loses its shear strength and liquefaction occurs. Liquefaction is 
the main reason for much of the destruction caused by the 
earthquake. Large and continuous deformations that occur in soil 
mass can lead to instability and ultimately great damage to these 
structureslMaurer, Green et al. 2014). This phenomenon can cause 
gable roofs instability, loss of bearing capacity of structures 
foundation, and damage to the structures caused by non-uniform 
subsidence. 

After the wide liquefaction that occurred in the devastating 
earthquake of Japan and Alaska in 1964, geotechnical engineers 
paid attention to this phenomenon. In the last 40 years, a significant 
progress has been done in terms of understanding liquefaction 
mechanism and the factors affecting it(Rahman, Siddiqua et al. 
2015). In the early years, more attention was focused on 
investigating this phenomenon in clean sands so that it was thought 
that liquefaction is just related to sands and coarse-grained soils are 
not capable to produce an additional pressure of pore water, which is 
the main reason of liquefaction (Seed and Idriss 1982). But as time 
passed and new earthquake occurred and following the observation 
of this phenomenon in coarse and fine-grained soils, many 
researchers aimed to examine the factors affecting the liquefaction 
of the soil (Idriss and Boulanger 2006). 

Iran is located on seismic Alpine-Himalayan region and this 
increased the occurrence of earthquakes as one of the most 
damaging natural disasters in the country. Major earthquakes such 
as Manjil earthquake (1990), Ardebil (1996), South Khorasan 


(1997), Bam (2003) and Boroujerd (2006) confirms this issue; 
therefore, the need is felt for a comprehensive review of earthquake 
phenomenon and its effects in the country. Therefore, in this study, 
Babol City (as shown in Figure 1) was evaluated in terms of 
geological and seismic features, using results of soil mechanic 
studies and eventually liquefaction risk zoning map of the city was 
provided(Choobbasti, Farrokhzad et al. 2015). 



Figure 1. Position of the Babol City in the Mazandaran Province 


2. METHODS 

Less than one century has passed from the first use of the word 
“liquefaction” by researchers. In fact, this is known as one of the 
most complex issues in geotechnical engineering. The Close 
connection of this phenomenon with earthquake engineering has 
attracted more attention from scientists. Hazen in 1918 used the 
word “liquefaction” for the first time in order to explain Calaveras 
dam break. Studies on liquefaction phenomenon were seriously 
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started by researches after two earthquakes of Niigata in Japan and 
Alaska in America in 1964 (Ceyhun and Hilmi). 

Among the most important studies in this field, I can refer to 
Casgrande, Castro, Seed and Idriss, Youd, Iwasaki et al.’s research 
works. In this study, according to Seed et al.’s method, the safety 
factor is determined against liquefaction in accordance with Eq. 1 
(Seed, Tokimatsu et al. 1984). 

CRR R Cyclic Resistance Ratio ( 1 ) 

5 CSR L Cyclic Stress Ratio 

Cyclic shear stress ratio is determined by Eq. 2. 

CSR = — = 0/65 — (2) 

t t d 

S °"o 

Where V is the average equivalent uniform cyclic shear stress 
caused by the earthquake, which is assumed to be 0.65 of the 
maximum induced stress; <2 max is the peak horizontal acceleration 
at ground surface generated by the earthquake, g is the acceleration 
of gravity, <T 0 and cr' 0 are the total and effective overloaded 

stresses, respectively, and T^ is a stress reduction coefficient (Seed 

and Idriss 1971). f d had been suggested by the national center for 
earthquake engineering research workshop (NCEER) in 1996 


(Youd, Idriss et al. 2001): 



r d =1.0-0.00765Z 

for Z <9.15 m 

(3) 

r d =1.174-0.0267Z 

for 9.15 m <Z <23 m 

(4) 

r d =0.744-0.008Z 

for 23 m <Z < 30 in 

(5) 

r d =0.50 

for Z )30 m 

(6) 


The first step in determining the cyclic shear resistance ratio 
(CRR) of the soil is the corrected number of SPT (N) to be 
achieved (iV^. 


(N ) -N C C C C C 

I iv D60 iv (7) 

Where (AM,„is the corrected number of SPT value, N is the 
standard penetration resistance measured at the site, C N is the 
correction factor N is considered as overhead effective stress, 

m 

C E is hammer efficiency, C B is the correction factor for borehole 
diameter, C R is the correction factor for bar length and C s is the 

correction factor for samples with or without cover. The NCEER 
workshop proposed the Equation 8 for fine correction coefficients 
(Youd, Idriss et al. 2001). 

( JV .)«C=“ + ^( JV .)«, < 8 > 

Where a and (3 are coefficients that are obtained based on the 
percentage of fines (F200). The cyclic shear resistance ratio (CRR) 
of soil is determined to obtain the value of (Nl)60cs and using the 
Seed et al.’s graphs with attention to the fine-grained soil (Youd, 
Idriss et al. 2001). In order to determine the safety factor against 
liquefaction in earthquakes with different magnitudes, the scale 
magnification correction coefficient of the earthquake must be used. 
In this study, the suggested coefficients of Poud & Idris are used. 


Poud and Idris suggested that MSF coefficients are calculated 
according to Equation 9. 

FS = (CRR 7.5 / CSR )MSF (9) 

Where MSF is the magnitude correction factor. 

3. DIFFERENT PROCEDURES FOR EVALUATING 
FINE-GRAINED SOILS LIQUEFACTION 

Natural sand deposits may include some silt seam and bands which 
can play an important role on the excess pore pressures being 
developed; therefore, this issue sometimes can be a critical factor in 
increasing the potential of liquefaction. Sites consists of silt and silty 
clays with low plasticity properties are widely available around the 
world, including Iran. Thus, it is very important to know the 
properties of silts encountered in field or used in laboratory testing. 
If the silt layer has a high plasticity, it can increase the resistance 
against liquefaction while low plasticity silts may alleviate the 
potential risk of liquefaction. 

According to Chinese criteria, sandy soils with low liquefaction 
properties have more trend to liquefaction at low shear strain and 
soils with high plastic properties tend to reduce the resistance and 
hardness at the large strain. Seed et al. provided new criteria in 2003 
for fine-grained soil liquefaction assessment (Seed, Cetin et al. 
2003). 

These criteria consider all the properties of fine-grained soil and are 
widely used in its engineering projects. This procedure is described 
in Figure 2. As it can be seen, if the soil is in A or B area, there is a 
need for more investigations including laboratory tests or confidence 
coefficient in SPT method. In order to calculate the liquefaction 
severity in this area, 3 methods of Haeri & Yasrebi, Iwasaki and 
Sunmez and Gukejlo were used. The procedure of Iwasaki includes 
liquefaction severity from ground level to a depth of 20 meters. 



LL (Liquid Limit) 

Figure 2. Recommendations Regarding of “Liquefiable” Soil Type 
Assessment 


In this way, by calculating the liquefaction potential at any location 
using equation 10, the effects of liquefaction in the study can be 
concluded (Iwasaki, Tokida et al. 1982): 


r 20 

/ L =j o FW (z)dz 

(10) 

fl-FS FS <1 


F=\ 

(ii) 

[0 FS )1 


W’(z) = 10-0.5z 

(12) 
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Where / is the liquefaction potential index. F is the function of the 

safety factor, FS is the safety factor against liquefaction at depth z 
and w [z ) i s the depth function. Iwasaki proposed four categories 

for assessing the severity of liquefaction. Liquefaction severity is 
very low for / =0- low for/ L (5 , high for 5(1 L (land very high 

for 15(7 L . 

Haeri and Yasrebi changed the safety factor in Iwasaki method, in 
order to remove defects in this approach in order to reduce the 
confidence coefficient effect of different layers against liquefaction 
in final potential function of surface occurrence. Also, the depth 
function was changed so that if the liquefaction layer depth 
increased, the effect of liquefaction in that layer would reduce non- 
linearly on the surface. Following equation is suggested to calculate 
the potential for liquefaction at surface(Iwasaki, Tokida et al. 1982). 


• 20 20-Z 
'■ 5 Z 


(1 -F)dz 


(13) 


Where P L the liquefaction potential index at surface, Z is the 
depth in meters, F is the safety factor against liquefaction at 
depth Z. 

Analyzes have been conducted by Haeri and Yasrebi show 
that P = 1 is the appropriate boundary for separating the 
two modes of occurrence and non-occurrence of liquefaction 
effects. Thus, the P L < lvalues indicate non-occurrence at 
surface and Pj > 1 are indicative of the occurrence of the 

surface effects of liquefaction(Haeri 1999). 

Chen and Zhang presented the possibility of soil 
liquefaction by PL equation. This index varies from zero to 
one which is associated with FL (Chen and Juang 2000). 



Sunmez and Gukejlo designated the LS liquefaction severity 
index instead of the liquefaction potential index proposed by 
Iwasaki. Sunmez and Gukejlo described the construction of 
liquefaction severity classified according to the following 
formula: 

L 5 =j o 2 V i ( z )( 1 0-0.5z)* (15) 

The PL (z) value is determined for Eq. 14 for each depth. 


L s =P L (z)( 10z 


z\|2o f 2 ° dP L (z ) 

4 1° dz 


(10z 



Classification of liquefaction severity index which was suggested by 
Sunmez and Gukejlo is given in Table 1 with a description of the 
liquefaction potential(Sonmez and Gokceoglu 2005). 


Table 1. Classification of the severity of liquefaction based on LS 
(Sunmez and Gukejlo) 


Severity of liquefaction 

Liquefaction Index (Ls) 

Very high 

85 <Ls< 100 

High 

65 < Ls < 85 

Medium 

35 < Ls < 65 

Low 

15 < Ls < 35 

Very little 

0<Ls< 15 


Non-liquefaction(FL>1.411) 


Ls=0 


4. EVALUATION OF THE GEOLOGICAL, SEISMICITY 
AND GEOTECHNICAL STUDIES AT BABOL CITY 

Babol, a city of Mazandaran province in the northern part of Iran, is 
considered as the study area in this research. The city is located 
approximately 20 km south of the Caspian Sea on the west bank of 
the river Babolrood and receives abundant annual rainfall. 

Babol area is located on alluvial deposits that belong to quaternary 
geology resulting from torrential rivers. The bedrock in the area is 
made of Conglomerate, sandstone. In general, the soil texture in this 
area consists of sand, clay and silt compounds. The percentage of 
these compounds depends on the distance from the sea. Thus, the 
closer to the sea, the higher percentage of sand compounds and 
contrariwise: the more the distance from the sea, the higher 
percentage of silt compounds. 

A large number of faults have been identified along the North East - 
South West direction intersecting with faults along the North West - 
South East in Babol region, which leads to smashing bedrock in this 
area. The acceleration of plan in this project is considered 0.35 g 
according to Babol location in a region with a relatively high risk. 
However, there is a considerable risk due to the Micro-faults in 
Quaternary sediments and smashing bedrock. The field of study is 
considered as a class (IV) according to the classification table of the 
new edition of the regulations [4]. 

According to the seismic investigations that have previously been 
conducted in the area of the Mazandaran province, the assessment of 
seismic hazard in the Babol city on the analysis of earthquake 
magnitude (M) equals to 6.4 on the Richter scale. This is obtained 
by the maximum horizontal acceleration (Amax) equal to 0.32g. 
This number roughly corresponded with the number achieved by the 
ministry of Housing and Urban Development in a 500 year-return 
period event for Babolsar area. The maximum horizontal 
acceleration for a return period of 500 years is estimated to be 
between 0.2g to 0.3g in Babol area. 

In a study that was conducted by American researchers to 
demonstrate the effect of site, it was proposed to the U.S. regulation 
draft, Martin and Dobry showed that for the accelerations less than 
0.4g, acceleration over alluvial is more than the acceleration of the 
stone, and to obtain the acceleration over alluvial it should be 1.2 
times the acceleration of the stone. From what has been discussed 
above, which were obtained from the latest studies and research 
works, the acceleration number (Amax) is obtained 0.35g in this 
research. Also, the earthquake magnitude has been considered 7.5 
on Richter scale for the calculations of liquefaction. Laboratory and 
geological studies imply the following results about Babol city soil: 

1- Geological formations of plan scope generally belong to the 
Pleistocene (the present era). 

2- Formations of the present era which has covered all over the 
studied fields have three origins of marshy-fluvial, and the flood 
and most of the three categories of sediments surround the 
sediments originating from the fenny and marshy areas. 

3- The dominant soil type is clay and sandy silt. 

4- Sand particles with a diameter of up to 0.2 mm, are the diameter 
of the dominant particles in sandy layers of the area. 

5- According to the USCS classification, the dominant soil type is 
the type of CL and ML (clay and silt with low pasty properties). 

6- Dry specific weight of the soil for the soil's area is about 1.1 
g/cm3 to 1.5 g/cm3. 

7- Wet specific gravity of the soils for the proposed area is located 
in the range of 1.6 g/cm3 to 2.1 g/cm3. 

8- Underground water depth varies from 0.2 to a maximum of 6.5 
meters in height in the more southern areas. 

Due to the city map, according to Figure 3 and considering the 
dispersion of boreholes location in the city, for appropriate 
conclusion and the best assessment of the area liquefaction, it was 
divided into zones of 1000 to 1,000 meters, in which 5 zones were 
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obtained in the horizontal direction and 7 zones in the vertical 
direction. Overall, the city was divided into 35 zones which can be 
seen in Figure 3. Overall, I will have 22 zones of important areas of 
the city for severity of liquefaction assessment that is usable for 
zoning map. 


I « I » I <= | o | . | 



Figure 3. Zoning map of the Babol city for liquefaction assessment 
and the boreholes position 

5. RESULTS AND DISCUSSION 


Table 2. Physical characteristics and soil testing of Borehole number 


33 


Borehole No. 3 Depth of groundwater level: 1.5 m 

Date of test: 91/05/18 

C(gr/cm2) 

-s 

W (%) 

D 50 (mm) 

P200 

Percent fine 

Twet 

(g/cm3) 

yd(g/cm3) 

Nspt 

Soil Type 

(USCS) 

Depth(m) 



27.63 

0.072 

64.5 

1.69 

1.32 

8 

ML 

2 



23.61 

0.121 

51.8 

1.69 

1.37 

8 

ML 

4 



28.45 

0.071 

65.0 

1.74 

1.35 

10 

ML 

6 

0.5 

32 

26.54 

0.148 

45.7 

1.82 

1.44 

14 

SM 

8 



24.66 

0.122 

49.1 

1.77 

1.42 

12 

SM 

10 



26.77 

0.122 

52.8 

1.64 

1.29 

7 

ML 

12 



22.65 

0.192 

40.9 

1.77 

1.44 

12 

SM 

14 



25.17 

0.116 

53.3 

1.8 

1.44 

13 

ML 

16 



24.37 

0.189 

38.2 

1.79 

1.44 

13 

SM 

18 



22.87 

0.198 

33.0 

1.77 

1.44 

12 

SM 

20 


In the fourth and fifth columns of Table 2, the special weight of dry 
and wet soil is provided. The parameter of wet special weight is 
used in calculation of effective and total tension. P 200 is the 
percentage of passing particles through the No. 200 sieve. This 
parameter is used in order to consider the amendments related to the 
assessed soil fine grained measure to determine the cyclic resistance. 
D 50 presents the particles diameter of soil that 50 percent of grains is 
small than it or in other words D 50 shows the average size of soil 
grain in millimeters. 


The physical properties and boreholes soil test were adjusted in the Table 3 shows the sam P le out P ut of written program for liquefaction 

table for easy access to the boreholes data. Table 2 also represents a potential assessment. This analysis was performed for all boreholes 

sample of these tables from No. 33 boreholes data. The results of and overa11 results related t0 the a 11 boreholes are visible in Table 4. 

liquefaction occurrence severity for each zone were calculated In this program, some soil parameters including special weight of 

separately with three methods of Iwasaki, Haeri and Yasrebi and sod and tbe percentage of fine-grained are known as a program 

Sunmez and Gukejlo, which were adjusted in Table 3. Finally, by input parameters. 


using a computer program written in this regard in Matlab Program, 
calculations were performed to evaluate the liquefaction hazard. 

Table 3. Sample output from a program written for liquefaction assessment 


Depth(m) 

2; 

•a 

F200 

Ol 

Ol 

3 

CSR 

CN 

(N1 )60 

(Nl)60c 

Soil Type 
(USCS) 

Depth(m) 

2 

8 

63.5 

3.28 

2.8 

0.98 

0.26 

1.7 

10.2 

17.24 

ML 

2 

4 

8 

52 

6.7 

3.2 

0.9 

0.35 

1.5 

9.81 

16.77 

ML 

4 

6 

10 

65 

10.2 

5.7 

0.95 

0.39 

1.32 

11.3 

18.46 

ML 

6 

8 

14 

45 

13.8 

7.2 

0.93 

0.4 

1.16 

15.4 

23.5 

SM 

8 

10 

12 

49 

17 

8.9 

0.9 

0.4 

1 

12.0 

19.48 

SM 

10 

12 

7 

52 

20 

10 

0.8 

0.39 

0.9 

6.93 

13.32 

ML 

12 

14 

12 

2152 

23.2 

11.7 

0.8 

0.38 

0.92 

11.0 

18.29 

SM 

14 

16 

12 

38 

27.8 

13.3 

0.75 

0.35 

0.87 

11.2 

18.51 

ML 

16 

18 

13 

42 

31.6 

14.9 

0.62 

0.34 

0.82 

10.63 

17.77 

SM 

18 

20 

12 

- 

33.9 

13.3 

0.63 

0.31 

0.78 

9.25 

15.91 

SM 

20 




X 

n 

Sunmez 

Index 







CRR 

Safety 

Factor 

Iwasaki 

Index 

o> 

3 . 

o' 

Q- 

X 

CD 

P 

Q- 

0> 

X 







0.11 

0.44 

10.7 

2.92 

0.09 

18.3 







0.1 

2 

0.0 

0.0 

0.4 

0.6 







0.12 

0.22 

10.2 

1.32 

0.99 

14.8 







0.16 

0.41 

7.66 

0.91 

0.98 

12.7 







0.1 

0.2 

7.4 

0.5 

0.9 

10 







0.1 

0.2 

7 

0.2 

0.9 

8.9 







0.12 

0.34 

4.71 

0.24 

0.99 

6.95 







0.12 

2 

0.0 

0.0 

0.03 

0.18 
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0.12 

2 

0.0 0.0 

0.04 

0.11 

0.11 

2 

0.0 0.0 

0.04 

0.04 

Table 4. Summary of Results Indicators boreholes and liquefaction 



zone 



Zone 

Borehole 

Iwasaki 

Haeri 

Sunmez 


Number 

Index 

Index 

Index 

A5 

45 

8.9 

0.5 

18.15 


18 

44.05 

2.63 

55.54 

B3 

32 

31.12 

4.88 

48.44 

29 

57.17 

6.11 

76.62 


Ave 

44.11 

4.54 

60.20 


35 

40.81 

91.19 

89.03 

B4 

50 

34.33 

6.4 

76.30 

34 

40.85 

5.68 

69.85 


Ave 

38.66 

7.09 

78.39 

B5 

11 

20.96 

5.49 

56.96 

B6 

25 

33.71 

5.23 

58.20 

B7 

36 

44.19 

3.32 

64.80 


23 

58.81 

7.87 

84.84 

Cl 

20 

21.93 

1.39 

32.29 


Ave 

40.37 

4.63 

58.57 


4 

24.11 

2.17 

43.23 

C4 

16 

50.29 

4.80 

68.72 

39 

57.46 

7.54 

84.73 


Ave 

43.95 

4.84 

65.56 


30 

16.41 

0.97 

26.46 


40 

16.47 

0.97 

26.47 

D2 

13 

32.10 

4.37 

57.55 


19 

55.95 

4.05 

69.93 


Ave 

30.23 

2.59 

45.10 


41 

62.27 

8.21 

89.77 

D3 

48 

28.73 

7 

55.87 


Ave 

45.50 

7.60 

72.82 


37 

53.48 

6.23 

80.89 


2 

23 

2.61 

35 

D4 

10 

55.65 

5.89 

82.17 


49 

51.85 

8.75 

89.69 


Ave 

46 

5.87 

71.95 


28 

40.70 

3.95 

63 

D5 

9 

0 

0 

2.66 


Ave 

20.30 

1.98 

32.86 


1 

28.14 

3.35 

44.34 

D6 

47 

25.64 

3.35 

49.54 


Ave 

26.89 

3.35 

46.94 


21 

36.63 

6.09 

59.92 


43 

43.15 

5.80 

63.38 

C2 

31 

33.95 

4.38 

52.23 


5 

32.32 

4.53 

53.58 


Ave 

36.51 

5.20 

57.53 


44 

23.84 

6.54 

57.20 


3 

56.21 

8.41 

87.97 

C3 

33 

47.68 

6.22 

73.85 

46 

40.01 

5.47 

57.18 


14 

61.76 

9 

97.42 


Ave 

45.90 

7.14 

74.72 


38 

68 

7.4 

98.8 

C5 

42 

62.33 

9.12 

98.27 


Ave 

65.19 

8.26 

98.54 


6 

61.23 

8.1 

88.46 

C6 

24 

26.31 

4.69 

45 


Ave 

43.77 

6.4 

66.74 

Cl 

15 

52.2 

5.53 

72.21 

D1 

12 

46.34 

8 

80.45 


26 

52.51 

5.29 

73.99 

D7 

17 

43.72 

4.4 

59.36 


Ave 

48.11 

4.85 

66.68 


8 42.9 3.51 64.56 

r i - 

27_48.76_ 2M _65 

Information of all data in Babol city can be seen in Table 4. Each 
borehole was analyzed by three approaches and liquefaction indexed 
was obtained for each borehole. To calculate the liquefaction index 
of different zones in this table, average of boreholes result was 
calculated in the last line of each zone where there was in that zone 
and finally three liquefaction index has been obtained for each zone 
based on Sonmez, Haeri and Iwasaki methods. 

In this map which obtained from Iwasaki procedure which is 
showed in Figure 4, areas are categorized into three diverse parts as 
low liquefaction potential, high liquefaction potential and very high 
liquefaction potential that almost all zones in Babol city was 
considered as a very high severity of liquefaction and just one zone 
was recognized as high liquefaction potential. Also there are four 
different liquefaction severities in Iwasaki’s approach; I can just two 
status of liquefaction as non-liquefaction and very high liquefaction 
potential. It seems that liquefaction will occur in almost all area of 
Babol city, based on obtained map from Iwasaki method 
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Figure 4. Microzonation Map of liquefaction severity of Babol City 
according to Iwasaki 



In Haeri and Yasrebi's approach, there are just two conditions 
for occurring or non-occurring liquefaction and the map 
which has been obtained from this approach are categorized 
into two sections as shown in Figure 5. In this method, almost 
all areas are recognized as liquefiable areas. In fact, Haeri 
method presents the presence or absence of liquefaction at 
surface, numbers of greater than 1 indicates presence of 
liquefaction in the surface and numbers of less than 1 
indicates absence of liquefaction in the surface. This map 
provides information about liquefaction segments like 
Iwasaki procedure. In this map 22 zones have been evaluated 
in terms of determining liquefaction and non-liquefaction 
areas. Overall, liquefaction will occur in 21 zones and this 
issue indicates that liquefaction index is greater than 1 in 
these zones and one zone is recognized as non-liquefaction 
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segment that can represent that liquefaction index is lower 
than 1. 

In Figure 6 which is obtained by Sonmez & Gokceoglu procedure, 
almost all existing zones on the map have been identified as areas of 
liquefaction except for zone E2, E3 and A6. On this map, just one 
zone (C6) in the center of the map is known as liquefaction zone 
with high intensity. Other areas that were identified with high- 
intensity liquefaction in previous procedures, in this method have 
been determined as liquefaction areas with high and medium 
intensity. Moreover, southern and central area of Babol are 
recognized as medium and high liquefaction potential whereas these 
segments were identified as very high liquefaction potential in 
Iwasaki approach. 
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Figure 5. Microzonation map of occurrence or non-occurrence of 
liquefaction in Babol city according to Haeri and Yasrebi approach 

(1976) 
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Figure 6. Microzonation map of liquefaction severity of Babol city 
according to Sunmez and Gukejlo (2005) 

But non-liquefaction areas have almost been the same in all maps 
and almost all 3 zones map, A3 and E2, E3 are known as non¬ 
liquefaction. This issue indicates that to evaluate the non¬ 


liquefaction areas, all three methods have acted in the same way, 
and the non-liquefaction areas are the same for all three approaches. 
By comparing the obtained map from the three proposed methods, it 
can be concluded that Iwasaki and Haeri’s method have an 
acceptable consistency because there are many similarities in the 
obtained maps from these two methods. On the other hand, the 
obtained liquefaction map from Sonmez’s approach gives different 
results. I can see similar results in Heidari et aKHeidari, Amel-Sakhi 
et ah), Farokhzad et al(Farrokhzad, Choobbasti et al. 2012) and 
choobbasti et al(Choobbasti, Farrokhzad et al. 2015) based on 
NCEER2001 and artificial neural network but the number data was 
not enough for receiving suitable results. In addition, for increasing 
the accuracy of results, the average of boreholes result was 
calculated for each zone. Finally, in this paper three maps were 
obtained from divers procedures were compared for enhancing 
validation of results. With evaluation and comparison of obtained 
maps with different methods can be concluded that the central and 
southern parts of the Babol city are known as liquefaction areas in 
almost all three maps, but there are some differences about northern 
parts of Babol city. Evaluation of obtained map from Haeri and 
Iwasaki method shows that all northern parts of Babol city are 
known as liquefaction point with high intensity, while in the 
obtained map from Sonmez procedure, northern parts of Babol city 
are known as liquefaction point with moderate intensity. However, 
about non-liquefaction areas, obtained maps have a lot of 
consistency. 

6. CONCLUSIONS 

The results of 50 bore holes from field experiments were used for 
preparing microzonation map of the city of Babol. The city has been 
divided to 35 separate zones for preparing microzonation map that 
22 zones were used for zoning. The zoning maps were gained from 
the adjusted results in table 4, in which the liquefaction intensity 
zoning map of Babol city is provided according to Iwasaki’s 
procedure. Concerning the map was obtained by this approach, 21 
zones were evaluated with very high liquefaction intensity and 1 
zone with high liquefaction severity from the set of 22 searchable 
zones. Obtained map from Haeri & Yasrebi’s approach indicates 
liquefied and non-liquefied segments in Babol. From a total of 22 
evaluable zones, the results show the liquefaction incidence at the 
surface for 21 zones and only one zone is determined with the lack 
of liquefaction incidence at the surface. The zoning map of 
liquefaction occurrence intensity in Babol city based on Sewmez 
and Gokojlou’s method illustrates that 2 zones with low liquefaction 
incidence severity, 9 zones with moderate liquefaction incidence 
severity, 10 zones with high liquefaction incidence intensity and 
only 1 zone with very high liquefaction incidence intensity were 
evaluated as shown in Figure 6. The zoning maps were obtained 
from Iwasaki and Haeri’s approach can confirm each other, so that 1 
zone with low liquefaction intensity according to Iwasaki procedure 
is the same zone which is evaluated without liquefaction occurrence 
at surface in Haeri & Yasrebi’s approach. Due to the different 
divisions of Sewmez & Gokojlou’s method or in the 2 previous 
procedures, I cannot have an appropriate analogy for results 
correspondence, hence; it is reasonable to conclude that the maps 
which obtained from Iwasaki and Haeri’s approaches give us more 
valuable results and I can trust on these maps. In addition, the main 
criteria of Iwasaki and Haeri’s procedure are so similar and these 
methods have more compliance with soil in Iran. 
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